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In maize developing seeds, transfer cells are prominently located at the basal endosperm 
transfer layer (BETL). As the first filial cell layer, BETL is a gateway to sugars, nutrients and 
water from mother plant; and anchor of numerous functions such as sucrose turnover, 
auxin and cytokinin biosynthesis/accumulation, energy metabolism, defense response, 
and signaling between maternal and filial generations. Previous studies showed that basal 
developing endosperms of miniature^ (mn^) mutant seeds lacking the Mnl-encoded 
cell wall invertase II, are also deficient for hexose. Given the role of glucose as one 
of the key sugars in protein glycosylation and proper protein folding; we performed a 
comparative large scale glycoproteome profiling of total proteins of these two genotypes 
(mnl mutant vs. Mn^ wild type) using 2D gel electrophoresis and glycosylation/total 
protein staining, followed by image analysis. Protein identification was done by LC-MS/MS. 
A total of 413 spots were detected; from which, 113 spots matched between the two 
genotypes. Of these, 45 showed >20% decrease/increase in glycosylation level and 
were selected for protein identification. A large number of identified proteins showed 
decreased glycosylation levels in mn^ developing endosperms as compared to the Mn^. 
Functional classification of proteins, showed mainly of post-translational modification, 
protein turnover, chaperone activities, carbohydrate and amino acid biosynthesis/transport, 
and cell wall biosynthesis. These proteins and activities were related to endoplasmic 
reticulum (ER) stress and unfolded protein response (UPR) as a result of the low 
glycolsylation levels of the mutant proteins. Overall, these results provide for the first 
time a global glycoproteome profile of maize BETLenriched basal endosperm to better 
understand their role in seed development in maize. 
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INTRODUCTION 

The mnl mutation in maize is associated with a loss of 70% 
seed weight at maturity due to a loss-of-function mutation at the 
Mnl locus that codes for the endosperm-specific cell wall inver- 
tase 2 (INCW2). The INCW2 protein is entirely localized in basal 
endosperm transfer cell layer (BETL), a major if not the sole gate- 
way for the intake of sugars and nutrients from maternal cells into 
the filial tissue (Miller and Chourey, 1992; Cheng et al, 1996). 
Homozygous mnl mutant is not lethal presumably due to a resid- 
ual low level of the cell wall invertase (CWI) activity encoded by 
another homolog, Incwl (Chourey et al, 2006). A BETL in devel- 
oping maize seed is composed of two to three strata of highly 
specialized transfer cells (Davis et al., 1990; Kang et al., 2009), 
marked by a hallmark feature of labyrinth-like proliferation of 
cell wall called wall-in-growths (WIGs) in plants (Offler et al, 
2003; Vaughn et al, 2007; Ruan et al, 2010). The INCW2 pro- 
tein is immunolocalized to the Mnl WIGs, and the mnl mutant 
is associated with stunted WIG development due to the INCW2 



deficiency (Kang et al., 2009). Sucrose hydrolysis is clearly a major 
physiological function of the Mnl -encoded INCW2 in the BETL 
and there is evidence that metabolically released hexoses, not 
exogenous, are critical as a driving force in assimilate movement 
between maternal and filial cells (Cheng and Chourey, 1999). As 
expected from the INCW2-deficiency, the mnl basal endosperm 
shows greatly reduced levels of glucose, fructose and sorbitol; and 
increased levels of sucrose relative to the Mnl developing seeds 
(LeClere et al., 2010; Chourey et al, 2012). 

Glycoproteins (GPs) have one or more covalently attached gly- 
can (oligosaccharide, often with the presence of glucose) moieties. 
Nearly all N-linked GPs are believed to be secretory proteins, 
which are associated with numerous functions, including stress 
tolerance and cell-cell communication in development. Proteins 
entering the secretory pathway are marked by signal peptides; 
17% of all proteins in Arabidopsis have predicted signal pep- 
tides and, 33% of these have at least one trans-membrane domain 
associated with the endoplasmic reticulum (ER). Addition of 
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W-glycans to the proteins occurs as a co-translational modifica- 
tion at tlie ER and then transported through secretory pathway 
to the cell surface in a traffic process mediated mostly by the 
Golgi complex (Vitale and Boston, 2008; Liu and Howell, 2010). 
The ER is suggested to be a nursery where newly synthesized 
secretory proteins are properly decorated, modified, de-decorated 
(deglycosylated), finally properly folded and assembled prior to 
exit to their final destination (Vitale and Boston, 2008). Quality 
Control (QC) in this entire process is highly critical. Improperly 
folded proteins in the ER are causal to the so called ER-stress 
that is associated with unfolded protein response (UPR) — a sig- 
nal transduction pathway that is well studied in mammalian and 
yeast cells. Recent studies in plants suggest that QC, ER stress, 
and UPR play an important role in biotic/ abiotic stress and seed 
development (Liu and Howell, 2010). 

The most detailed studies of GPs in maize developing 
endosperm have been on the storage protein zein, stored in ER- 
derived protein bodies (Vitale and Boston, 2008; Arcalis et al, 
2010). However, there is no such information on the GPs in 
the BETL, a cell type not associated with storage functions in 
seed development. Here we report the first detailed profile of the 
GPs in the basal endosperm using LC-MS/MS approach, a pre- 
eminent tool for identification and quantitative characterization 
of GPs (Ruiz-May et al., 2012). Comparative GP profiles of the 
BETL enriched proteins of the mnl relative to the M«l revealed 
potential changes due to both the hexose deficiency. 

EXPERIMENTAL PROCEDURES 
PLANT MATERIALS AND CHEMICALS 

Immature maize (Zea mays L.) kernels wild type {Mnl) and 
miniature 1 mutant (mnl) in the W22 inbred line were har- 
vested at 12 days after pollination (DAP). All plants were grown 
in the field and were self- or -sib-pollinated. At the time of 
harvest, kernels were individually excised from the ear with 
a paring knife, taking care to include undamaged base (pedi- 
cel) of each kernel. Excised kernels were flash frozen in liquid 
nitrogen and stored at — 80°C until analysis. We used the basal 
1/3 end of the endosperm because it is enriched for the BETL 
cells, the sole site of the Mnl expression and also a major 
zone for sucrose turn-over reactions, as previously discussed 
(LeClere et al., 2008; Silva-Sanchez et al., 2013). All chemicals 
were purchased from Fisher Scientific Inc., USA unless otherwise 
stated. 

TOTAL PROTEIN EXTRACTION 

Soluble proteins were extracted from the one-third lower part 
of maize kernels according to a method reported by Hurkman 
and Tanaka (1986) with minor modifications. The frozen ker- 
nels were ground into fine powder in a pre-chilled mortar and 
pestle. A total of 3 mL of extraction buffer (0.1 M Tris HCl pH 
8.8, 10 mM EDTA, 0.2 M DTT, 0.9 M sucrose) and 3mL of sat- 
urated phenol were added to homogenize the sample. To extract 
proteins, samples were constantly agitated for 2 h at room temper- 
ature. The mixture was then centrifuged at 5000 x g for lOmin 
at 4°C. The top clear phenol phase was removed and collected 
in a fresh tube. The remaining pellet was re-extracted with 3 mL 
buffered phenol solution and centrifuged again to recover the 



top clear phenol phase and combined with the previous fraction. 
The supernatant was precipitated with five volumes of ice cold 
0.1 M ammonium acetate in 100% methanol overnight and cen- 
trifuged at 20000 X g for 20min at 4°C. The supernatant was 
discarded and the pellet was washed twice with 0.1 M ammo- 
nium acetate in 100% methanol, then washed twice with cold 
80% acetone and finally with 70% of ethanol. The pellet was 
solubUized in isoelectric focusing (lEF) buffer (7M urea, 2M 
thiourea, 4% (w/v) CHAPS and 40 mM DTT). The samples were 
treated with benzonase (Novagen, Gibbstown, NJ) for 30min 
and then centrifuged at 34,000 rpm for 30min at 15°C using an 
ultracentrifuge (Optima TLX ultracentrifuge, Beckman Coulter). 
Supernatant was collected and 50 |xL aliquots were kept at — 80°C 
until use. 

TWO DIMENSIONAL GEL ELECTROPHORESIS, PROTEIN STAINING AND 
IMAGING ANALYSIS 

Protein samples were quantified using an EZQ protein quantita- 
tion kit (Invitrogen, CA, USA). Protein extracts of 500 (xg each 
were mixed with 0.2% of ampholytes (pH 3-10) and loaded 
onto an 18 cm Immobilized pH gradient (IPG) strips (pH range 
3-10 NL) (GE healthcare, CA, USA), followed by overnight rehy- 
dration at room temperature. lEF was conducted on an Ettan 
IPGphor3 system (GE Healthcare, CA, USA) using the following 
conditions: 200 V for 30 min, then ramping to 500 V for 30 min, 
and finally to 10,000 V for 1 h. The voltage was held at 10,000 V 
until 85,000 Vh were reached in order to ensure complete sep- 
aration and focusing of the proteins. After lEF, the strips were 
equilibrated for 15 min in equilibration buffer (6 M urea, 75 mM 
Tris-HCl pH 8.8, 29.3% glycerol, 2% SDS, 0.002% bromophenol 
blue) containing 2% DTT, and for another 15 min in equilibra- 
tion buffer containing 2.5% iodoacetamide. The IPG strips were 
placed onto 18 cm 12.5% SDS gels (Jule Biotechnologies INC, 
Milford, CT, USA). Electrophoresis was run at 15 W for 5h. 
A total of three replicate 2 -DE experiments were conducted for 
each sample. 

After gel electrophoresis, staining with a glycoprotein-specific 
stain Pro-Q Emerald fluorescent dye was performed according 
to manufacturer's instructions (Invitrogen, CA, USA). Images 
were acquired using the Investigator ProPic unit equipped with 
a UV light-box (Genomics solutions, MI, USA). Then, the gels 
were destained and restained with a total protein dye Sypro 
Ruby following manufacturer's instructions (Invitrogen, CA, 
USA). Images were acquired with the Investigator ProPic unit. 
Spot detection, matching and quantification across the repli- 
cate gels of wild type and mutant samples, were done using 
a Progenesis Samespot Software (Non-linear Dynamics, CA, 
USA). The function of automatic spot detection and match- 
ing was used, followed by manual inspection. Protein experi- 
mental molecular weights were calculated using a CandyCane™ 
molecular mass standard (Invitrogen, CA, USA) separated in 
a separate lane on the 2D gels, and the experimental isoelec- 
tric points were determined based on IPG strip specifications. 
Normalized spot volumes were used to determine the quanti- 
tative changes of glycosylated/total protein ratios and total pro- 
tein levels. ANOVA test was used to determine the statistical 
significance. 
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PROTEIN IDENTIFICATION USING LIQUID CHROMATOGRAPHY 
TANDEM MS (LC-MS/MS) 

Selected 2D gel spots were excised using an Investigator ProPic 
robot (Genomics Solutions Inc., USA) and digested with trypsin 
as previously described (Sheffield et al, 2006); The lyophilized 
peptides were resuspended in 15 (xL of loading buffer (3% ace- 
tonitrile, 0.1% acetic acid, 0.01% trifluoroacetic acid) and loaded 
onto a C18 capillary trap cartridge (LC Packings, USA) and then 
separated on a 15 cm nanoflow analytical C18 column (PepMap 
75 |xm id, 3 |i,m, 100 A) at a flow rate of 300 nL/min on a nanoLC 
ultra ID plus system (ABsciex, USA). Solvent A composition was 
3% acetonitrile (ACN) v/v, 0.1% acetic acid v/v; whereas solvent 
B was 97% ACN v/v, 0.1% acetic acid v/v. Peptide separation 
was performed with a linear gradient from 3 to 40% of solvent 
B for 20 min, followed by an increasing to 90% of solvent B in 
5 min and hold for 5 min (Zhu et al., 2010). The eluted peptides 
were directly sprayed into an LTQ Orbitrap XL mass spectrome- 
ter (Thermo Scientific Inc., Bremen, Germany). MS2 spectra were 
acquired in a data-dependent mode. An Orbitrap fall MS scan 
(resolution: 3 x lO'*, mass range 400-1800 Da) was followed by 
10 MS2 scans in the ion trap, which were performed via col- 
lision induced dissociation on the top 10 most abundant ions. 
Isolation window for ion selection was 3 Da. Normalized colli- 
sion energy was set at 28%. Dynamic exclusion time was 20 s (Li 
et al., 2012a). The acquired mass spectra were searched against 
a UniprotZeci mays database (62,860 entries, Jan-04-2013) using 
Mascot 2.2 search engine (http://www.matrixscience.com) with 
the following parameters: tryptic peptides with 1 missed cleav- 
age site, mass tolerance of precursor ion of lOppm and MS/MS 
ion of 0.8 Da, fixed carbamidomethylation of cysteine, variable 



methionine oxidation, asparagine and glutamine deamination. 
Unambiguous identification was done using Scaffold software V 
3.0 (Proteome Software, OR, USA). 

BIOINFORMATICS ANALYSIS 

Functional classification of the proteins was done according to the 
clusters of orthologous groups for eukaryotic complete genomes 
(KOG) (ftp://ftp.ncbi.nih.gov/pub/COG/KOG/) (Tatusov et al, 
2003) and the prediction of subcellular localization was done 
using a Plant-mPLoc tool (http://www.csbio.sjtu.edu.cn/bioinf/ 
plant-multi/) (Chou and Shen, 2010).The identified proteins 
were submitted to an AT-Glycosite tool for glycosylation site anal- 
ysis (http://www.hiv.lanl.gov/content/sequence/GLYCOSITE/ 
glycosite.html) (Zhang et al., 2004). Prediction of the N-terminal 
presequences: chloroplast transit peptide (cTP), mitochondrial 
targeting peptide (mTP) and secretory pathway signal peptide 
(SP) was done using the TargetP 1.1 tool (http://virww.cbs.dtu. 
dk/services/TargetP/) (Emanuelsson et al., 2000). 

RESULTS AND DISCUSSION 

TWO DIMENSIONAL GEL ANALYSIS, PROTEIN IDENTIFICATION AND 
GLYCOPROTEIN (GP) DETECTION 

In order to determine the differences in glycosylation patterns 
of BETL between the wild type (WT) and the mutant (M), 
total protein extracts were prepared from the one-third lower 
part of maize kernels. For each sample, three biological repli- 
cates/gels were performed. The gels were stained with ProQ 
Emerald for glycoproteins, and then for total protein detection 
with Sypro Ruby stain. Figure 1 shows representative images of 
gels stained for potential glycoproteins and total proteins. A total 
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FIGURE 1 I Representative 2D-PAGE of BETL total protein extract 
stained with ProQ emerald for glycoprotein detection (A) and Sypro 
ruby for total protein detection (B). Tine spots along with the spot 



numbers selected for protein identification by LC-MS/MS are marked with 
a circle. Numbers on various spots in each gel refer to the #s shown in 
the 1st column, Table 1. 
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of 413 spots were detected in the four gels of the three replicates 
(Supplemental Figure 1) representing the wild type and mutant 
endosperm extracts; of these, 254 spots with p-values smaller than 
0.05 based on ANOVA test were selected for further analysis. 

Glycosylation of each protein spot was evaluated by their 
staining intensity with the GP stain that was normalized by 
the intensity of total protein stain (G/TP) for both WT and 
the mutant samples. Normalized glycosylation ratios are rep- 
resented as M(G/TP)/^T(G/TP) in Table 1. A total of 113 pro- 
tein spots matched between the wild type and mutant gels, are 
shown in Supplemental Table 1 . There are 45 spots (Figure 1) 
showing more than 20% increase or decrease in the normal- 
ized glycosylation level. Interestingly, only 4% of the 113 spots 
matched between WT and M showed increased GP ratios when 
comparing to the WT, while 36% showed decreased ratios 
(Supplemental Table 1). 

Figure 2 shows four representative spots with dynamic 
changes in normalized glycosylation levels (G/TP). The majority 
of proteins have decreased glycosylation ratio that is calculated 
by M(Gijp)/'WT(Qjjpy, as demonstrated by spots 357 and 436 
in Figure 2 with ratios of 0.69 and 0.71, respectively. Only four 
spots showed 20% increased glycosylation in the mutant com- 
pared to the wild type. The spot 547 in Figure 2 is one of those 
proteins with glycosylation levels in the mutant being more than 
two-fold of those in the wild type. This indicates the relation- 
ship between glycosylation and protein abundance is complex. 
Figure 2 also includes an example of a protein, spot 185, with the 
glycosylation ratio of 0.97 that represented no significant change 
in glycosylation (Supplemental Table 1) in the two genotypes. 

The 45 spots from the wild type and the mutant were excised 
and subjected to trypsin digestion and LC-MS/MS analysis. After 
database search using Mascot, the data were filtered in Scaffold 
software, and protein identification was considered valid if it had 
an overall 99% confidence and at least three unique peptides 
with 95% confidence. A total of 35 spots showed 69 compara- 
ble protein identifications between the WT and theM (Table 1). 
Some spots such as 272, 436, and 449 yielded multiple protein 
identifications due to the overlapping nature of proteins with 
similar molecular weights and isoelectric points (pis). Several 
proteins were also reported by Zhu et al. (2006) in a similar 
proteomic study of root elongation zone in maize. Conversely, 
the same proteins were identified from multiple spots, e.g., spots 
425, 454, and 461 as subtilisin-like protease, spots 357, 362, and 
371 as Rhicadhesin receptor, and spots 440 and 449 as ER lumi- 
nal binding protein. Differential post-translational modifications 
may account for these isoforms (Satoh et al., 2002). Alternative 
splicing is known to create multiple RNAs from a single gene 
that increase the proteome diversity. In maize, 19.2% of the total 
expressed genes are reported to be alternatively spliced (Barbazuk 
et al, 2008); thus various protein isoforms encoded by the same 
gene or its homologs may have alternative sizes and/or pis. 

Although nearly all the identified proteins (Table 1) exhib- 
ited close correlations between experimental and reported pis, 
some did show discrepancy between the experimental molecular 
weights and theoretical molecular weights (e.g., spots 376, 454, 
and 493). These variations can be attributable to several factors, 
such as the lack of full length cDNA sequences or misannotation 



for these proteins; some of the cDNA sequences may have differ- 
ent splicing or degraded products, and post-translational mod- 
ifications like glycosylation, which can also alter the molecular 
weights. Out of the 69 proteins presented in Table 1, 39 have been 
reported to be glycosylated. 

FUNCTIONAL CLASSIFICATION 

Functional classification of the identified proteins (Table 1) was 
performed according to the clusters of orthologous groups 
for eukaryotic complete genomes (KOG, Tatusov et al., 2003). 
Thirteen different categories were found: (I) Post-translational 
modification, protein turnover, chaperone functions; (II) carbo- 
hydrate metabolism and transport; (III) Amino acid metabolism 
and transport; (IV) Energy production and conversion; (V) 
Translation; (VI) RNA processing and modification; (VII) 
Nucleotide metabolism and transport; (VIII) Signal transduction; 
(IX) Coenzyme metabolism; (X) Secondary metabolites biosyn- 
thesis, transport and catabolism; (XI) Defense mechanisms; (XII) 
General functional prediction only; (XIII) No KOG related; as 
shown in Figure 3A. 

Category I (Post-translational modification, protein turnover, 
chaperone functions) was the most represented category that 
included 15 proteins such as ER luminal binding protein (spots 
440 and 449, Table 1), glutathione transferase (spot 362), heat 
shock proteins (spots 449 and 513), proteasome subunits (spots 
376 and 433), chaperone proteins (spot 449), vignain (spot 
440), subtilisin like protease (spots 454, 461, and 425), and WD 
repeat (spot 491). Interestingly, all proteins represented in the 
10 spots are glycosylated as reported previously in other organ- 
isms (Table 1, the last column for the references); however, all 
showed reduced normalized glycosylation ratio except for the heat 
shock protein 70 (spot 513) with a ratio of 1.19. These spots 
showed a total protein ratio, M{jp)/WT(jp), of 1.0 or slightly 
higher (Table 1), suggesting that the reduced glycosylation in the 
mutant was not associated with the reduced abundance of these 
proteins. 

Protein glycosylation is accomplished when the lipid 
(dolichol) -linked oligosaccharide transfers the carbohydrate- 
containing structure (Glc3Man9GLcNAc2) to nascent peptides 
in the ER (Howell, 2013). Several studies have shown that the 
glucose residues on the lipid linked oligosaccharide facilitate 
the in vitro transfer of the oligosaccharide to the protein 
although it is not an absolute requirement (Kornfeld and 
Kornfeld, 1985). Cells starved of glucose or in energy depriva- 
tion state fail to produce Glc3Man9GlcNAc2-P-P-Dol; instead, 
they accumulate Man9GlcNAc2 lipid linked species that were 
associated with the decrease in protein glycosylation in thy- 
roid slices (Spiro et al., 1983). The accumulation of defective 
saccharide-lipid linked carrier upon different stresses varies 
according to cell type. In the BETL region of the mnl mutant 
kernels, the low glucose content may lead to defective or 
low levels of the dolichol-linked oligosaccharide causing a 
reduction of the normal glycosylation process associated as a 
general trend of low glycosylation ratios when compared with 
the WT. 

Table 1 depicts proteins with a high total protein ratio when 
comparing to the wild type, including heat shock proteins (spots 
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FIGURE 2 I Examples of protein spots from BETL total protein extracts 
[spot 185 lM^G/TP)/WT^G/TP) = 0.97], spot 357 lMfG/TP)/VVT(a/m = 
0.74], spot 436 [M^GfTP)/WTfG/TP) = 0.69], and spot 547 IM^g/tP)/ 
WT(G/TP) = 2.35] showing quantitative differences in glycosylation (G) 
and total protein (TP) between WT and mutant after glycoprotein and 



Sypro Ruby staining. Tine spots tliat sliowed a 20% increase/decrease 
(357, 436, and 547) in tine glycosylation ratio were identified by mass 
spectrometry and labeled with the same numbers as in Table 1. Spots 
that did not showed significant changes in glycosylation ratios (185) are in 
Supplemental Table 1. 



449 and 513), glutathione transferase III (spot 362) and pro- 
teasome (spots 376 and 433). Proteomic studies in developing 
seeds of Brassica campestri L reported that 23.4% of the total pro- 
teins are involved in protein processing and degradation (Li et al., 
2012b). 

One of the main protein folding machineries in protein pro- 
cessing is the N-glycan dependent folding pathway (Howell, 
2013). After the Glc3Man9GlcNAc2 core is transferred to the 
nascent polypeptide by the oligosaccharide transferase (OST), the 
Glc a-1, 2, and a- 1,3 are removed by the glucosidase I and II, 
then calreticulin/calnexin along with protein disulfide isomerase 
(PDI) help in the proper folding of the protein. Folded proteins 
are released by glucosidase which cleaves off the last glucose a- 
1,3 and are finally secreted. Proteins that are not properly folded 
are sensed by UDP-glucose:glycoprotein glucosyltransferase, and 
they can be reglucosylated to reenter the calreticulin/calnexin 
mediated folding cycle; where N-glycosylation plays an impor- 
tant role in the endoplasmic reticulum quality control (ERQC). 



We hypothesized that ER stress was induced in the BETL region 
due to the observed low protein N-glycosylation, and since 
N-glycans are recognized along several steps in the folding pro- 
cess of proteins, low N-glycosylation may lead to accumulation 
of unfolded/misfolded proteins that triggers the UPR response 
to mitigate the ER stress. These observations are consistent with 
Howell (2013) that UPR in plants is associated with upregulation 
of genes involved in protein folding and endoplasmic reticulum 
associated degradation (ERAD). 

We observed four Hspl and Hsp70 proteins (spots 449 and 
513); some with a high protein ratio in the mnl mutant kernels 
when compared with the WT; these proteins are reported to be 
involved in protein targeting to the mitochondria and chloro- 
plasts (Kriechbaumer etal, 2012), along with ER luminal binding 
proteins (spots 440 and 449), proposed to be involved in pro- 
tein body assembly within the ER (Li et al, 2012b), supporting 
the hypothesis of ER stress and UPR process as a response of the 
defective glycosylation process occurring in the BETL region. 
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FIGURE 3 I Functional classification (A), and subcellular localization prediction (B), of the 69 proteins identified by LC-MS/MS in total protein 
extracts. 



In addition to UPR response, the degradation of proteins via 
the ubiquitination-mediated pathway is an important mecha- 
nism in plant growth and development, and is responsible for the 
degradation of abnormal peptides and short-lived cellular regu- 
lators; controlling many processes that allow rapid response to 
intracellular signals (Sadanandom et al., 2012).This pathway is 
also related the ERAD by sensing glycoproteins that fail to fold 
into a native state. The alpha mannosidase interfere the reglucosy- 
lation cycles and removes the a- 1,6 linked mannose from the core 
oligosaccharide linked to a glycoprotein. The mannose trimmed 
glycoprotein is recognized by a series of luminal proteins such as 
Lectin OS9 and HRD3 and then exported to the cytosol to be 
degraded via ubiquitinilation mediated pathway (Howell, 2013). 
The observation of proteasome subunit alpha and beta (spots 433 
and 376 respectively) in a high total protein ratio in the mnl 
mutant compared to the WT kernels suggest that mnl kernels are 
undergoing ER stress due to low concentration of glucose in the 
BETL region. 

A correlation between ER stress and programed cell death 
(PCD) is also described in plants (Zuppini et al., 2004). PCD 
occurs when internal contents of a cell are engulfed in the vac- 
uole, leading to the tonoplast rupture and release of vacuolar 
hydrolytic enzymes. The tonoplast rupture is mediated by vac- 
uolar processing enzymes that have caspase-like activity (Howell, 
2013). Along with the observed proteins related to ER stress, UPR, 
and ERAD processes, we found that the subtilisin-like protease 
(spots 425, 454, and 461) also have a reported caspase-like activity 
(Vartapetian et al., 2011). This protein showed a low normalized 
glycosylation level and no change in the total protein ratio in the 
mnl mutant compared to the WT, suggesting that ER stress might 
regulate PCD through under-glycosylation in the mnl mutant 
BETL region. Overall, the reduced glycosylation levels of many of 
the above proteins may also affect diverse metabolic and cellular 
functions in developing seeds. 

Carbohydrate cataboiism and transport, the second most rep- 
resented group in the functional classification (Figure 3A), is 
composed of eight proteins in the glycolysis pathway, pentose 



pathway, fructose and mannose metabolism and starch synthesis 
(spots 272, 403 with two proteins, 425, 436, 590, 626, and 649). 
The enzymes that participate in glycolysis pathway showed low 
ratios of glycosylation in the mnl mutant when compared to WT 
and no change in the total protein ratio (Table 1). Several gly- 
colytic enzymes among other protein/enzymes are reported to be 
modified by O-GlcNAc in various systems (Love and Hanover, 
2005). The enzymes that participate in glycolysis are of special 
interest, due to the greatly reduced levels of hexoses resulting from 
the INCW2- deficiency in the BELT region of the mutant. 

One of the proteins involved in carbohydrate cataboiism 
is glucose-6-phosphate-isomerase (spot 425). It catalyzes the 
conversion of glucose-6-phosphate to fructose 6-phosphate, an 
important step in the hexoamine biosynthetic pathway. We 
observed that the glycosylation ratio of glucose-6-phosphate- 
isomerase was reduced (~0.76) and the total protein ratio was 
nearly the same (~1.16), suggesting that its enzymatic function 
and activity may be modulated through glycosylation. Cells can 
sense the level of a product, UDP-GlcNAc, in the hexoamine 
biosynthetic pathway as an indication of the level of glucose, the 
precursor that is used in the same pathway (Love and Hanover, 
2005). Because glucose-6-phosphate-isomerase is an essential 
enzyme in the hexoamine biosynthetic pathway, we hypothesize 
that it plays an important role in maintaining the homeostasis in 
the cells in terms of nutrient levels in the mnl kernel. Further 
studies are needed in the characterization of this particular pro- 
tein and its role in regulatory pathways. 

Triosephosphate isomerase (spot 649) that catalyzes the inter- 
conversion of glyceraldehyde-3-phosphate to dihydroxyacetone 
phosphate was found in cytosolic and plastid isoforms with 
the same trend in glycosylation and total protein abundance as 
the glycolytic enzymes, in which glycosylation level goes down 
while the total protein level doesn't change significantly. For 
fructose and mannose metabolism, fructokinase-2 (spot 272) 
was found with a low glycosylation ratio but a high ratio in 
total protein (~1.5); while diphosphate-fructose-6-phosphate 1- 
phosphotransferas (spot 590) was found with a high ratio of 
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glycosylation (~1.2) and no significant change in the total pro- 
tein ratio. The 6-phosphogluconate dehydrogenase (spot 436), a 
decarboxylating enzyme, participates in the pentose phosphate 
pathway (PPP) and was found with low ratios for both glyco- 
sylation and total protein. Significantly, Spielbauer et al. (2013), 
reported a subset of chloroplast-localized PPP enzymes to be 
not only present in endosperm amyloplasts but are also critical 
for starch and oil biosynthesis in developing seeds. The UTP- 
glucose-1 -phosphate uridylyltransferase (spot 403) is responsible 
for synthesis and pyrophosphorolysis of UDP-glucose, a key pre- 
cursor of carbohydrate formation (including sucrose, cellulose, 
starch, glycogen and P-glucan biosynthesis). This protein showed 
one of the highest ratios in glycosylation but had a reduced ratio 
in total protein level. UGPase is highly abundant in cell wall mem- 
brane in barley and it is speculated to have a role in providing 
UDP-glucose for biosynthesis of cell wall components, including 
P-glucans and cellulose (Eimert et al., 1996). 

Amino acid metabolism and transport was the third most 
important category in functional classification. There were six 
proteins involved in amino acid metabolism; 2-isopropylmalate 
synthase, alanine aminotransferase, chorismate synthase, ser- 
ine hydroxymethyltransferase, transaminase and glutamine syn- 
thetase (spots 495, 436, 187, 590, and 547 for last two proteins). 
The first three enzymes showed a low glycosylation ratio but 2- 
isopropylmalate synthase and chorismate synthase had a high 
protein abundance ratio (Table 1). The last three enzymes had 
high glycosylation ratios but did not change in protein levels, 
which suggests that glycosylation may play a role in the regulation 
of these enzymes. In rice seed development proteomics studies, 
Deng et al. (2013) observed a correlation between the accumu- 
lation of proteins involved in protein synthesis and turnover, 
protein folding and amino acid metabolism, in early stages of 
rice seed development, indicating an active turnover of proteins 
during early development stages; as observed in our work. 

SUBCELLULAR LOCATION OF THE IDENTIFIED PROTEINS 

When observing the subcellular location prediction of the com- 
mon proteins between the wild type and the mutant (Figure 3B), 
the chloroplast was the most represented with a 33% contribu- 
tion of the total proteins identified, followed by the cell wall 
proteins (15%) and cytoplasm (13%). The chloroplast protein 
group consists of a very wide spectrum of different proteins that 
are involved in energy production and conversion, amino acid 
metabolism and transport, nucleotide metabolism and transport, 
carbohydrate metabolism and transport, and post-translational 
modification-protein turnover-chaperone functions. 

All cell wall associated proteins showed low glycosylation level 
in the mnl mutant compare to the WT, as expected. The cell wall 
group was composed of 10 proteins; three subtilisin-like proteases 
that have post-translational modification, protein turnover, and 
chaperone functions (spots 425, 454, and 461) did not show any 
significant change in their total protein ratio. The other seven 
proteins did not belong to any KOG related group and were 
germinin-like proteins (GLPs), rhicadhesin receptor, and alpha- 
JV-arabinofiiranosidase A (spots 357, 361, 362, 371, 449, 451, and 
498), It has been reported that about 40% of all known GLPs 
are cell wall proteins; and all GLPs contain N-terminal secretory 



sequences that relate them strongly to the cell wall/extracellular 
matrix targeting (Breen and Bellgard, 2010). GLPs have been 
studied as defense molecules in different plant species, conditions 
and diseases due to their high resistance to protease, heat, SDS 
and extreme pH. Several reports identified GLPs in diverse envi- 
ronmental conditions, such as salt, aluminum, drought stresses, 
fungal pathogen attacks, bacteria and virus infections (Breen and 
Bellgard, 2010). A special class of GLPs that is almost exclusively 
found in cereals is characterized by an oxalate oxidase (OXO) 
activity; meanwhile, a barley GLP reported by Bernier and Berna 
( 200 1 ) presents a different ADP glucose enzyme activity that leads 
to the hypothesis of its involvement in the control of metabolic 
flow toward starch, cell wall polysaccharides, glycoproteins and 
glycolipids in plants (Rodriguez-Lopez et al., 2001). In addi- 
tion, Gucciardo et al. (2007) isolated a pea GLP that showed 
similarity in its N-terminus with a rhicadhesin recepetor. The 
GLP mRNA was expressed in non-conventional locations for rhi- 
cadhesin receptors such as in nodules and the expanding cells 
adjacent to the nodule meristem and in the nodule epidermis. 
This protein is demonstrated to have a superoxide dismutase 
activity and show resistance to high temperature among other 
stresses (Breen and Bellgard, 2010). AH GLP proteins showed an 
increased protein ratio in the mnl mutant kernel, suggesting an 
upregulation due to the stress caused by the glucose-deficiency in 
BETL region. 

Eight proteins of cytosolic location (spots 362, 403, 418, 
425, 440, 480, 491, and 547) are described in Table 1 and 
Figure 3B were found in the experiment. These proteins have a 
variety of functions such as amino acid metabolism and trans- 
port; carbohydrate metabolism and transport; post-translational 
modification, protein turnover, chaperone functions; secondary 
metabolites biosynthesis; transport and catabolism; and trans- 
lation. Glutathione transferase III (A), peroxidase 2, and a like- 
WD-repeat protein (spots 362, 418, and 491; respectively) showed 
low normalized glycosylation level but an increase in the total 
protein ratio. These proteins are related to oxidative stress. In a 
previous study (Silva-Sanchez et al, 2013), proteins related to 
gluthatione metabolism were up-regulated, in the mnl kernels, 
suggesting a potential redox regulation by glutathione. Elongation 
factor 2, alcohol dehydrogenase class-3 and glucose-6-phosphate 
isomerase (spots 425, 440, and 480; respectively) showed a low 
normalized glycosylation level and nearly no change in their total 
protein ratio. These observations suggest that their regulation 
may be mediated through their glycosylation. A protein similar 
to UTP-glucose-1 -phosphate uridylyltransferase (UGPUT) from 
Zea mays (B6T4R3, spot 403) showed a high normalized gly- 
cosylation ratio (~1.48) but a low total protein ratio (~0.76). 
This protein plays an important role in the biosynthesis of cell 
wall polysaccharides as well as in the synthesis of the carbohy- 
drate moiety for glycolipid and glycoproteins (Kleczkowski et al, 
2004). Glutamine synthetase (GS) (spot 547) showed the highest 
level of normalized glycosylation (~2.35) but no change in total 
protein ratio (~1.01). Cytosolic GS participate in seed develop- 
ment by assimilating ammonium in sink tissues probably through 
the asparagine catabolism as a possible source of ammonium 
(Bernard and Habash, 2009). The high glycosylation levels of 
UGPUT and GS suggest a complex mode of post-translational 
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regulation in the plant. It is interesting to add that glycosylation 
of cytosolic proteins has been reported previously. Haltiwanger 
et al. (1992) described that several cytoplasmic and nuclear pro- 
teins undergo O-glycosylation, vifhich is highly dynamic in nature 
and is known to response to various stimuli in the media, in a sim- 
ilar way to phosphorylation responses. These processes are thus 
far well known in yeast or animal models (Spiro, 2002; Love and 
Hanover, 2005) although there is some research done in plants 
(Thornton et al., 1999; Love and Hanover, 2005 and references 
there in). 

/V-GLYCOSYLATION PREDICTION AND SECRETORY PATHWAY 

For all of the proteins present in Table 1, we predicted the num- 
ber of probable JV-glycosylation sites based on the N-X-S, N-X-T, 
or NX[ST]Z motifs, where X can be any amino acid residues 
except proline using the N-glycosite tool (Zhang et al., 2004). A 
total of 48 proteins showed at least one predicted N-glycosylation 
site in their sequence. Those proteins that did not show any 
predicted N-glycosylation sites potentially have other types of 
unusual or rare glycosylation or are not studied/reported as gly- 
coproteins yet. Although some proteins that showed at least one 
predicted N-glycosylation site, they may have been reported with 
other types of glycosylation. For example, the proteasome alpha 
subunit (spot 433) has been reported to be O-glycosylated in 
mouse (Overath et al., 2012); elongation factor 2 (spot 480) is an 
O-glycosylated protein in human cervical cancer cells (Solorzano 
et al, 2012), triosaphosphate isomerase (Spots 626 and 649) (Love 
and Hanover, 2005); and protein Z showed O-glycosylation in 
human plasma (Spiro, 2002). O-glycosylation occurs in a great 
variety of proteins that are involved in key nuclear and cytoplas- 
mic activities. Similar to phosphorylation, O-GlcNAc modifica- 
tion occurs on Ser or Thr residues. O-GlcNAc is interchangeable 
with phosphorylation. These two modifications have similar roles 
in biological functions as well. 0-GlcNAc modification is asso- 
ciated to a variety of regulatory activities such as carbohydrate 
metabolism, signaling, transcription and translation, and stress 
response (Love and Hanover, 2005). 

In order to investigate whether the proteins are related to 
the secretory pathway and have the potential of being glycosy- 
lated or not, the identified proteins from the 45 spots in Table 1 
were subjected to signal peptide prediction using the Target? tool 
(Emanuelsson et al, 2000). According to the prediction, 17 pro- 
teins showed the highest scores for SP and have signal peptides. 
Examples of these proteins included rhicadhesin receptor (spots 
357, 362, and 371), ER luminal binding protein (spots 440 and 
449), GLP's (spots 361, 451, and 498), subtUisin-like protease 
(spots 425, 454, and 461). Many of these proteins are cell wall 
proteins or membrane proteins and are known to be glycosylated 
(Denecke et al, 1991; Bykova et al, 2006; Koseki et al, 2006; 
Gucciardo et al, 2007). For proteins with relative low scores for 
signal peptide prediction, eight proteins were assigned to mito- 
chondria while 10 were located in the chloroplasts. The majority 
of these proteins have not been reported as glycoproteins. For 
the proteins that were not assigned neither to mitochondria or 
chloroplast; the low score on the SP prediction could indicate that 
these proteins have no classical secretion signal (Radhamony and 
Theg, 2006), and may follow so-called unconventional secretion 



based on alternative models of secretion independent of ER and 
Golgi route (Zhang and Schekman, 2013). 

CONCLUSIONS 

The study of glycosylated proteins in the total protein extracts 
of BETL from wild type and the mnl mutant revealed the defi- 
ciency of the mutant in the glycosylation process. The identified 
proteins that participate actively in post-translational modifica- 
tions, protein turnover and chaperone functions may indicate 
the ER stress and UPR in the mutant due to the defective gly- 
cosylation processes. Carbohydrate metabolism related proteins 
showed in general low glycosylation ratios, suggesting that their 
enzymatic functions and activities may be regulated by the gly- 
cosylation levels through the hexosamine pathway. Low glycosy- 
lation ratios were also found in cell wall proteins of the mutant, 
which represent most of the plant defense proteins and structural 
protein-crosslink. More studies are needed to fully understand the 
mechanism underlying protein glycosylation, sugar sensing and 
plant survival. 
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